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Verfahren zur Verbesserung der Tiefendiskrimiiiierung optisch-abbildender Systeme 

Die Erfindung betrifft ein Verfahren zur Tiefendiskriminierung optisch abbildender Systeme. 
Es ist insbesondere in der Lichtmikroskopie anwendbar um eine Erhohung der Bildqualitat 
5 bei der Untersuchung dreidimensional ausgedehnter Objekte zu erreichen. 

In der klassischen Lichtmikroskopie besteht bei der Untersuchung dreidimensional 
ausgedehnter Objekte, <Lh. von Objekten, deren Ausdehnung entlang der optischen Achse 
groJJer ist als die Tiefenscharfe der verwendeten Objektive, das Problem, dass das scharfe 
10 Bild aus dem Fokusbereich mit aufierfokalen Bildanteilen, welche demzufolge unscharf 
abgebildet werden, uberlagert ist. 

Um dieses Problem zu tiberwinden ist die Nutzung der konfokalen Abbildung bekannt, bei 
der mittels eines sogenannten Pinholes das aus dem nicht im Fokus liegenden Bereich 
stammende Licht ausgeblendet wird und somit nicht zur Abbildung beitragt. Auf diese Weise 
15 entsteht ein sogenannter optischer Schnitt Diese konfokale punktweise Abbildung erfordert 
das Abrastern des Objektes in der Bildebene um ein Bild zu erhalten. Diese Abtastung kann 
entweder mittels Scannern in Laser Scanning Mikroskopen oder mittels Nipkow-Scheiben 
erfolgen. 

Durch Aufhahme mehrerer optischer Schnittbilder in verschiedenen Fokuslagen kann ein „z- 
20 Stapel" gewonnen werden, welcher eine dreidimensionale Darstellung des Objektes 
ermoglicht 

Ein anderer Weg zur Erzeugung optischer Schnitte ist die Anwendung der strulcturierten 
Beleuchtung. Diese wurde erstmals von Meir Ben-Levy und Eyal Pelec in WO' 97/6509 
angegeben. Verbesserungen imd Erweiterungen dieser Methode wurden von Tony Wilson et 
al. in WO 98/45745 und von Volker Gerstner et al. in WO 02/12945 beschrieben. Auf die 
Offenbarung dieser drei Schriften wird hiermit als Referenz verwiesen. 
In WO 97/6509 wird das Objekt mit einer periodischen Struktur (Sinus- oder Rechteckgitter) 
beleuchtet und mittels einer Kamera ein Bild des Objektes aufgenommen, digitalisiert und in 
30 einem Speicher abgespeichert AnschlieBend wird die periodische Struktur innerhalb der 

Bildebene so verschoben, dass die Phasenlage der Struktur geandert wird, und wieder ein Bild 
aufgenommen und abgespeichert. Dieser Vorgang (Verschiebung, Bildaufeahme, 
Abspeicherung) kann mehrfach wiederholt werden. Anschlieflend werden die Bilder 
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miteinander verrechnet um ein Schnittbild zu erzeugen. Als mathematischer Ansatz wird 
dabei eine Fourierentwicklung angegeben, was zu einem komplizierten Formelapparat fuhrt. 
In WO 98/45745 wird eine einfachere Formel fur die Schnittbilder angegeben, welche sich 
durch Vereinfachung des Ansatzes aus WO 97/6509 fur den Falls jeweils gleicher 
5 Phasenverschiebungen zwischen den einzelnen Aufeahmen herleiten l&sst. 

Die Realisierung der fur die Anwendung der angegebenen Verfabren notwendigen 
physikalischen Randbedingungen erweist sich in der Praxis als sehr schwierig. So fuhrt z. B. 
die Schwankung der Lampenhelligkeit zwischen den verschiedenen Aufeahmen zu 
streifenformigen Artefakten in den generierten Schnittbildern. Bei Fluoreszenzobjekten treten 
10 zusatzliche Probleme durch die zeitabhSngige Ausbleichung der Fluoreszenzfarbstoffe auf, 
welche gleichfalls zu Fehlern ftihrt Auch die notwendige Konstanz der einzelnen 
Phasenverschiebungsschritte l&sst sich in der Praxis nicht einhalten. 
In WO 02/12945 wird daher vorgeschlagen, mittels Auskopplung eines Teils des zur 
Beleuchtung des Objekts dienenden Lichtes, Registrierung der Intensitat und anschliefiende 
15 Normierung der Einzelaufhahmen den Einfluss der zeitlich unterschiedlichen 
Lampenhelligkeit auszugleichen: Fur die Berixcksichtigung ungleicher 
Phasenverschiebungsschritte wird ein Gleichungssystem (dort GL 22) angegeben. Zum 
Ausgleich der Ausbleichung wird vorgeschlagen, statt der mindestens notwendigen drei 
Aufeahmen je Schnittbild 6 Aufeahmen in der Reihenfolge 1-2-3-3-2- Izu 
20 registrieren, zwei Schnittbilder (aus 1-2-3 und 3 - 2 - 1) zu berechnen und von diesen 
beiden den Mittelwert zu bestimmen. 

Die Realisierung dieser Vorschlage erfordert erheblichen instrumentellen Aufwand. 
AuBerdem verlangert die Aufeahme zusatzlicher Bilder die notwendige Aufeahmezeit und 
erhoht damit auch die Belastung der Probe durch Beleuchtung mit dem 
Fluoreszenzanregungslicht 

Die Erfindung stellt sich die Aufgabe, die Nachteile des Standes der Technik zu vermeiden 
und ein verbessertes Verfahren zur Bestimmung optischer Schnittbilder anzugebeiL 
Diese Aufgabe wird durch ein Verfahren zur Verbesserung der Tiefendiskriminierung optisch 
abbildender Systeme gemaiJ dem Hauptanspruch gelost Vorteilhafte Ausgestaltungen sind in 
3 0 den Unteranspruchen angegeben. 

Die Erfindung wird im folgenden anhand der Figuren naher beschrieben, 
Es zeigen: 

Fig. 1 den prinzipiellen Aufbau eines Mikroskops mit strulcturierter Beleuchtung 
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Fig. 2 eine schematische Darstellung der erfindungsgem&Ben Bandpassfilter in der 
Fourierebene 

Fig. 3 schematisch die Lage der Abbildungen der periodischen Struktur 

Fig. 4 einen Ablaufplan fur die Bestirnmung der Korrekturwerte durch lineare Optimierung. 

5 

In Fig. 1 ist das vereinfachte optische Schema der strukturierten Belenchtung am Beispiel der 
Durchlichtbeleuchtung dargestellt. Eingezeichnet ist der bildgebende Strahlengang 
(Lukenstrahlengang). Eine 1-dimensionaleperiodische Struktur (Transmissionsgitter) (3), 
welche sich in einer Scharfeebene der dargestellten optischen Anordnung befindet, wird durch 

10 eine Lichtquelle (1) und nachfolgender Kollektoroptik (2) beleuchtet Dem Gitter folgt in 
Lichtrichtung eine planparallele Glasplatte (4). Der Winkel der planparallelen Platte zur 
optischen Achse kann definiert eingestellt werden. Durch die nachfolgende 

^ beleuchtungsseitigen Optiken (5, 6) (Kondensor) wird die Struktur in die Praparatebene (7) 
H abgebildet. Von dem Praparat ausgehendes Licht wird wiederum iiber ein Linsenpaar (8, 9) 

15 (Objektiv und Tubuslinse) in die darauffolgende SchSrfenebene (10) abgebildet, in welcher z. 
B. eine CCD-Matrix einer digitalen Kamera angeordnet sein kann. Dabei dient die 
planparallele Glasplatte (4) durch definierte Verkippung zur Verschiebung der Abbildung der 
Gitterstruktur (3) auf das in der Praparatebene (7) befindliche Objekt. Vorzugsweise dient im 
Fall der Auflicht-Fluoreszenz-Beobachtung das Objektiv (8) gleichzeitig als Kondensor. 

20 Wie in WO 02/12945 beschrieben wurde werden jetzt mit mindestens 3 Stellungen der 

Glasplatte (4) mittels der digitalen Kamera jeweils eine Registrierung der sich ergebenden 
Helligkeitsverteilung im Objekt durchgefuhrt. 

Dabei kann die Helligkeitsverteilung Ij (x, y) fur den Fall eines Sinus/Kosinus-Gitters 
vereinfacht durch die Formel 



angegeben werden, wobei if=0 ... N-l die i-te Phasenposition des projizierten Gitters und N die 
Anzahl der Registrierungen ist, m(x, y) die Modulationstiefe des Objekts (und damit die 
30 gesuchte Bildinformation an der Stelle x, y) ist und die §\ die Phasenwerte sind. 

Diese Gleichung enthalt die drei unbekannten GroBen I 0 , m und (|)o. Somit kann durch 
mindestens drei Messungen bei gezielt variiertem §\ (i = 1,2,3) die Bestimmung dieser GrSBen 
erfolgen. Die Losung kann durch einen lcleinste-Quadrate-Ansatz aus den Messungen 




J, (x, y) = I 0 (x, y) - (1+ mix, y) - cos(& (x 9 y) + ^ (x, y)) 



(1) 
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gewonnen werden. Dazu bietet es sich an, Gleichung (1) in einer kompakteren Form 
darzustellen und die Cosinus-Funktion mittels Additionstheorem umzuschreiben: 

K*> y,</>i) = a 0 y) + a x (*, y) • f x ) + a 2 (x, y) • f 2 (4 ) (2) 
mit 

/ 1 (^) = cos^ i 
/,(^) = sin^ 

ao(*>J>) = /o(*>J>) 

«i (*> y) = (*> y) • ™(*> ■ c o s ( x > y) 

<* 2 (*> y) = ~ 7 o ^) • rn{x, y) • sin 0 O 7) 



Die Funktionen fi und f 2 h&ngen somit nur von den im Prinzip firei zu wahlenden 
Phasenverschiebungen <|>j ab. Fur den Fall nichtsinusformiger, aber periodischer Struktur kann 
die Helligkeitsverteilung Ii (x, y) auch durch Reihenentwicklung approximiert werden, das 
Prinzip der Rechnung bleibt dabei im Grande gleich. Die kleinste-Quadrate-Losung lautet in 
Matrixschreibweise 



M-a=b 



mit 



(3) 



M = 



(4) 



mit N = Anzahl Messungen (hier: Phasenschritte) und 



c = 



'a 0 (*,j>) N 



(5) 
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sowie 



E/<*,*#i) 

f t [l(x,y^ l )-fM)] 



V N 



(6) 



5 Die Losung dieser Matrixgleichung wird mathematisch durch invertieren der Matrix M 
gewonnen. 
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= AT' -b 



(7) 



Somit konnen aus den Messungen, dargestellt durch den Vektor b, und den 
Phasenverschiebungen des Gitters die gesuchten unbekannten Grofien ermittelt werden. 
Insbesondere ergibt sich die Modulationstiefe am Ort (x, y) zu 



15 
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(8) 



Da mit Defokussierung des Objektes der Gitterkontrast und damit zwangslaufig die 
Modulationstiefe im Objekt abnimmt, steilt Gleichung (8) eine Bestimmungsgleichung zum 
Erhalt eines tiefendiskriiiiinierten Bildes (optischer Schnitt) dax. Gleichungen (2) und (3) 
liefern einen universalen LSsungsweg, der die Anpassung der Phasenlagen §\ an die jeweils 
experimentellen Gegebenheiten erlaubt. So ist etwa die Wahl der Zahl der Schritte und die 
Positionen im Objektraum (Ort der abgebildeten Gitterlinien) frei wahlbar. 



Jedes optische System ubertragt gemaB seiner optischen Transferfunktion Raumfrequenzen 
bis zu einem Grenzwert, der durch die numerische Apertur NA der Optik und der 
25 Lichtwellenlange X bestimmt ist: 



AtzNA 
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Dies hat zur Folge, dass von einem in den Objektraum abgebildeten Gitter, welches kein 
reines Sinus-Gitter ist, die hoheren Hannonischen der Gitterperiode bis zu eben dieser 
Grenzfrequenz kmax tibertragen werden und sich in dem Fourier-transformierten Bild des 
5 Objektes als lokale Maxima wiederfinden. Auch in einem nach dem oben angefuhrten 

Schema ermittelten Schnittbild sind diese Komponenten des Gitters enthalten und resultieren 
in Streifenartefakten. Diese Artefakte sind dadurch eliminierbar, dass die Harmonischen des 
projizierten Gitters durch Fourier-Transformation im Raumfrequenzspektrum lokalisiert und 
durch gezieltes Setzen von Filtem (Bandpass-Filter, Bandstops) entfernt werden. Erne 
10 anschlielJende Rucktransformation fuhrt dann zu artefaktfreie Bildern. Diese Methode ist 
schematisch in Fig. 2 dargestellt. In der Fourierebene (11) sind mit (12) die hoheren 
Harmonischen der Gitterstruktur (3) bezeichnet, (13) bezeichnet den Nullpunkt und entspricht 
dem Gleichanteil der Beleuchtung, d.h. einer gleichmaBigen, nicht strukturierten Beleuchtung. 
Durch Einbringen den Harmonischen (12) entsprechender Bandpass-Filter in den 
15 Strahlengang tragen diese Bereichen nicht zur Abbildung bei und werden damit nicht von der 
digitalen Kamera registriert. 

Bei fluoreszierenden Objekten tritt ein weiterer Grund zum Auftreten von Streifenartefakten 
hervor, namlich das Fluoreszenz-Bleichen (photobleaching, fluorescence fading). Hierbei 
kommt es durch die Abbildung der Struktur (z.B. des Gitters) in das Objekt zu Regionen 
20 unterschiedlich starker Ausbleichung der Fluoreszenz und somit schliefllich zu 

Streifenartefakten in jedem in der Folge aufgenommenem Bild. Die Frequenzanteile der 
Artefakte im Spektrum eines solchen Bildes sind besonders stark bei der Grundfrequenz der 
projizierten Struktur. Durch Entfernen der Raumfrequenzanteile der projizierten Struktur im 
Fourier-transformierten Bild durch entsprechende Bandpass-Filter und anschlieBender 
Riicktransformation kann auch hier prinzipiell ein artefaktfireies Bild erzeugt werden. 

Die Verhinderung von diesen bei der Fluoreszenzbeleuchtung auftretenden Artefakten kann 
auch durch die Bestimmung und Ausnutzung der Bleicheigenschaften von fluoreszierenden 
Objekten bei der sequentiellen Aufhahme von Bildern desselben Objekts erreicht werden. 
30 Hierbei ergibt sich pro Bildaufeahme durch die Beleuchtung eine Bleichung des Objelcts an 
den beleuchteten Stellen. Die Starke des Bleichens hangt neben der Beleuchtungsintensitat 
von den photochemischen und photophysikalischen Eigenschaften des fluoreszierenden 
Objektes ab. Die Bestimmung der Fluoreszenzstarke des Objelctes durch Integration der 
Intensitaten tiber alle oder ausgewahlte Objektorte (x, y) ermoglicht die Reskalierung der 
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Bildsequenz nach erfolgter Aufhahme der Bilder. Dazu werden entweder uber analoger 
Kamera und Analog-Digital- Wandier oder einer CCD-Kamera in einem Computer 
gespeicherten Bildinfonnationen durch paarweise Ermittlung von Quotienten der 
Fluoreszenzstarken der Bilder die Intensitaten an den Bildorten (x,y) auf einen Maximalwert 
normiert. Somit lasst sich der Fluoreszenzverlust durch Fluoreszenzbleichen beriicksichtigen. 
Dazu konnen nach Fig. 3 Regions of Interest (ROT) definiert werden, die bei beiden 
BMaufhahmen gleichermaflen beleuchtet werden und in Folge fluoreszieren, das MaB der 
Abnahme der Fluoreszenz in diesen Bereichen ist ein MaB ftir die Bleich-Rate und kann damit 
fur die Normierung der verschiedenen Registrierungen der Helligkeitsverteilung 
herangezogen werden. Die Definition dieser ROIs kann dabei sowohl vom Benutzer iiber 
entsprechende Eingabemittel als auch automatisch (aus der jeweiligen Lage der abgebildeten 
Struktur) erfolgen. In Fig. 3 stellt (14) die Verhaltnisse bei einer ersten Registrierung der 
Helligkeitsverteilung auf dem Objekt dar, (15) die Verhaltnisse bei einer zweiten 
Registrierung mit veranderter Phasenlage der abgebildeten Struktur (16). Die Regions of 
Interest (17) werden vorteilhafterweise so definiert, dass in beiden Registrierungen 
gleichmaBige Beleuchtung erfolgt. 

Das erfindungsgemaBe Verfahren erlaubt die Bestimmung von optischen Schnitfbildern 
praktisch ftir beliebig gewahlte Phasenverschiebungen §\ . Zur Verhinderung von Fehlern sind 
die Phasenverschiebungen §\ mit hoher Genauigkeit zu bestrmmen. Dieses lasst sich in 
einfacher Weise mit dem im folgenden beschriebenen Kalibrierungsverfahren sichern, bei 
dem die Einflusse des gesamten Systems beriicksichtigt werden. 

Dazu wird zunachst ein Referenzbild (Objekt mit xiberlagertem Gitter) aufgenonnnen und 
registriert. Daraufhin erfolgt eine geringe Phasenverschiebung des Gitters auf eine zunachst 
unbekannte Phasenlage des Gitters (es soli diese „Verschiebung pro Ansteuerungssignal" erst 
kalibriert werden) und erneute Bildaufhahme. Die erhaltenen zwei Bilder werden nun 
rechnerisch verglichen. Dies kann ohne Einschrankung der Allgemeinheit durch 
Differenzbildung, Summierung oder jede andere Art erfolgen, die auf eine Merit-Funktion 
fiihrt. Resultat ist zunachst im Allgemeinen ein Streifen-behaftetes Bild. Die Schritte: 
Verschiebung des Gitters urn einen ldeinen Betrag, Bildaufhahme, Vergleich mit dem 
Referenzbild werden so lange wiederholt, bis die Merit-Funktion den Extremwert erreicht, 
also z.B. ein Differenzbild nur noch aus Untergrundrauschen besteht oder ein Summenbild 
maximale Werte erreicht. Ist letzteres erfUllt, ist die Kalibration abgeschlossen und ein so 
erhaltener Ansteuerwert kann in einem Speichermedium zum spateren Zweck abrufbar 
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gehalten werden. Dieser Wert entspricht damit genau der Verschiebung des Gitters urn eine 
voile Periode. Alternativ ist es auch moglich, ein Summeribild auf gleichmaBige 
Helligkeitsverteilung (d.h. Verschwinden der Streifenstruktur) zu bewerten. In diesem Fall 
entspricht dieser Wert dann der Verschiebung des Gitters um eine halbe Periode. 
Vorzugsweise wird diese Prozedur mit einem Spiegel als Objekt durchgefiihrt. Der 
beschriebene Vorgang kann sowohl manuell als auch automatisiert erfolgen. 

Im folgenden wird an Hand der Fig. 4 ein alternatives Verfahren zur Bestimmung der 
Einfliisse von Lampenhelligkeitsschwankungen, Ausbleichen des Objektes bei 
Fluoreszenzheleuchtung und iucht-sinusformigem Gitter beschrieben. 

Dazu wird die zeitliche Veranderung der Beleuchtung durch intensitatsproportionale skalare 
Faktoren 0 f > 0 beschrieben werden, welche bei jeder aufeinanderfolgenden Aufhahme 

berucksichtigt werden (die Abhangigkeit der GroBen gi, Oi, usw. vom Ort (x,y) wird im 
1 5 Folgenden zur Vereinfachung nicht mehr explizit angegeben): 

Si=Ofli (9) 
Die modellierte Beobachtung g/ ergibt sich nun aus dem Produkt mit dem idealen, 
konsistenten Phasenbild o,. Um vollstandig ein gestortes System zu beschreiben, braucht man 
nur N-l Faktoren 0 t zu bestimmen. 

20 

Das Ausbleichen des Objekts bei Fluoreszenzbeleuchtung wirkt in Abhangigkeit von der 
einwirkenden Strahlungsintensit&t liber die Zeit. Diese ist wegen des projizierten Gitters 
drtlich verschieden und erfordert deshalb eine Behandlung, die die Gitterfunktion 
berucksichtigt. Eine einfache Beschreibung wird moglich, wenn man das Ausbleichen als eine 
lineare Funktion der Beleuchtungsintensitat annimmt: 

k, - (1- rfo M k H jQ^k yi iZl mit k- 0 =1 (10) 

Der aus ortsvarianten Falctoren bestehende Vektor 0 < k, ^1 soil abhangig vom Verlauf der 
• Bestrahlungsintensitat und einem wahlbaren Mafi 1 > d ^0 fur das Ausbleichen die 
Abschw&chung beschreiben wenn man ihn komponentenweise mit den jeweils noch nicht 
30 ausgeblichenen Beobachtungen 0/ multiphziert: 

g, =o,k, (11) 

Aus (10) geht hervor, dass immer jede Beobachtung mit den Falctoren aller vorhergehenden 
Phasenaufhahmen multipliziert wird, diese sind bereits z.T. ausgeblichen und xiberlappen sich 
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. mit der momentanen Phase. Altemativ kann man k* auch exponentiell d.h. entsprechend der 
Abklingkurve eines bestimmten Fluorochromes beschreiben. 

Nichtsinusffirmige Gitterverlaufe verursachen Streifenartefakte insbesondere fur ganzzahlige 
5 Vielfache der Gitterfrequenz. Unter genauer Kenntnis von Frequenz, Phasenlage und 

Kurvenform des Gitters bei den einzelnen Aufiiahmen wird es moglich einen sinusformigen 
Verlauf zu erzeugen. Ein Korrekturvektor l Kor r, ahnlich wie unter (10), hier jedoch entlang der 
Gitterperiodizitat, soil unter (12) betrachtet werden. 

1^=7^7- (12) 

10 Die Nachbildung der observierten Gitterfunktion ist mit faiuer gekennzeichnet Ein ideal 

sinusformiger Intensitatsverlauf wird mit fsmus zum Ausdruck gebracht. Die Konstante s im 
Nenner soil die Division durch sehr kleine Werte von f enter iibemehmen. Die Nachbildung der 
observierten Gitterfunktion kann auf verschiedenen Wegen geschehen. Eine einfache 
Moglichkeit einer Approximation bietet sich z.B. durch die Verwendung einer Fourierreihe 
15 fur die Trapezfunktion an. 

r / i_x 1 2 Asin(2fa*-&) . _ . N 
2 bn " (2z —1) 

Mit b kann man die Flankensteilheit der Trapezubergange variieren, x ist die Ortskoordinate 
und M die Anzahl der benutzten Fourierkoeffizienten. Alternativ lassen sich auch genauere 
Modelle verwenden. Diese erfordern die Faltung der wahren Gitterfunktion mit einer lateralen 
20 Komponente der PSF des Mikroskops. Von entscheidender Bedeutung ist eine phasen- und 
frequenzsynchrone Synthetisierung des Korrekturvektors l^ orr . Ist dies gewahrleistet, so kann 
ahnlich wie unter (11) mit einer einfachen komponentenweisen Multiplikation ein korrigiertes 
Ergebnis fur die Observation erzielt werden. 

fZlKorr ~ § i^lKorr (14) 

25 Die korrigierten Phasenbilder erhalt man nun aus: 

/,=o,==-^^,mit 0,>Onk,>O (15) 

Durch Einsetzen von (15) in Gleichung (3) bzw. (8) erhalt man eine allgemeine Formel zur 
Rekonstruktion, welche auBerdem die Parameter 9 i , d und b enthalt 

30 Die Losung des gestellten Problems kann beispielsweise durch Umformung in eine 

Extremalaufgabe erfolgen. Mit den Hilfsmitteln numerischer Optimierung kann man die 
Parameter 0 t 9 d und b unter Zuhilfenahme einer Merit-Funlction bestimmen. 

9 
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M(0 l9 d;b) « a 0 \F{5} 0 \ 2 + a I |F{2}.|* +a 2 \F{a} 2<D \ 2 + ... + a n \F{a} nto \ 2 -> min (16) 

F{ a } ist die Fourier-Transformation einer oder mehrerer Komponenten des Resultatvektors 5 
aus Gleichung (5). Dabei sind neben der Fourier-Transformation auch andere 
5 Funktionaltransformationen geeignet 

Ziel ist es 0 ( , d und b so zu variieren, dass (16) minimiert wird. Dafiir bietet sich eine 
Vielzahl von Methoden an die keinesfalls eingeschrankt zu Gradienten- oder 
Liniensuchmethoden zu sehen sind. Jeder Transformationskoeffizient ist auflerdem gewichtet 
10 mit einem Koeffizienten Damit kann der Algorithmus an verschiedene Signal-Rausch 
AbstSnde oder bevorzugte Frequenzen angepasst werden. Ein diesbezuglich vorteilhafter 
Wert den Gleichanteil zu wichten ergibt sich mit folgender Beziehung: 
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lF{5} 1 l 2 + lF{a} 2 l 2 +... + |F{5} n | 



a, oc r giii ui^ii ■ - ■ r (17) 



Die allgemeine Funktionsweise des hier beschriebenen Korrekturapparates ist dem 
schematischen Ablau^lan in Fig. 4 zu entnehmen. Die Verwendung von Projektionen ist 
vorteilhaft um effizient die Parameter 0 t , d und b zu bestimmen. Die Parameterermittlung ist 

aber nicht auf diesen eindimensionalen Fall eingeschrankt, sondern funktioniert in der Ebene 
20 in analoger Weise. 

Es werden bei fester Fokuslage von der Kamera mindestens drei Helligkeitsverteilungen (18), 
(19), (20), ggf. auch mehr (21), mit verschiedenen Phasenlagen i=l, 2, 3, .., N der auf das 
Objekt abgebildeten Gitterstruktur registriert. Mittels Projektion werden diese 2- 
dimensionalen Helligkeitsverteilimgen jeweils in eine 1-dimensionale Verteilung (22), (23), 
(24), und ggf. (25) uberfiihrt. Mit Anfangswerten fur die Parameter 0 t , d und b werden nach 

Formel (15) korrigierte Phasenbilder berechnet und aus diesen mittels des Ansatzes aus 
Formel (3) ein erstes Bild des optischen Schnitts errechnet (26). Die Anfangswerte der 
Parameter 0 t , d und b konnen dabei vom Bediener eingegeben oder auch automatisch 

gesch&tzt werden. F(ir diese ScMtzung konnen bekannten Tools der Minimierung / 
30 Maximierung einer Merit-Funktion benutzt werden (z. B. Simplex, Marquardt, o.a.). Aus 

diesem Ergebnis wird nun in (27) die Merit-Funktion nach Formel (16) bestimmt und in (28) 
auf das Erreichen eines Minimums getestet. Ist das Minimum noch nicht erreicht wird eine 
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neue Schatzung der Parameter 9 t , d und b durchgefuhrt (29) und mit (26) .fortgesetzt. Fiir den 
Fall des Erreichens des Minimums der Merit-Funktion, d.h. weitere Variation der Parameter 
0 { , d und b lassen ihren Wert wieder ansteigen, sind die optimalen Parameter bestimmt und 

in (30) wird das optische Schnittbild mittels der Formeln (15), (3) und (8) bestimmt und in 
(3 1) ausgegeben. AnschlieBend kann die Fokuslage verandert werden und nach dem 
dargestellten Verfahren ein weiteres optisches Schnittbild bestimmt werden usw. . Durch 
dieses Abrastern des Objektes in z-Richtung lasst sich ein z-Stapel von Schnittbildem und 
damit eine dreidimensionale Darstellxing des Objekts generieren. 

Die Realisierung der Erfindung ist nicht an die dargestellten Ausfuhrungsbeispiele gebunden, 
durch fachmannische Weiterentwickiungen oder andere Ansatze zur Beriicksichtigung der 
Helligkeitsschwankiihgen der Lichtquelle oder des Ausbleichen des Objektes bei der 
Fluoreszenzbeleuchtung wird der Bereich der Erfindung nicht verlassen. 
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1. Verfahren zur Verbesserung der Tiefendiskriminierung optischer Systeme bei dem 

a) ein Objekt mit einer periodischen Struktur beleuchtet wird, 

b) eine Registrierung der so entsteheriden Helligkeitsverteilung erfolgt, 

c) die Phasenlage der periodischen Struktur verschoben wird, 

d) die Schritte a) bis c) sooft wiederholt werden bis mindestens drei 
Helligkeitsverteilungen registriert wurden, 

e) die registrierten Helligkeitsverteilungen miteinander verrechnet werden urn eine 
Objekthelligkeitsverteilung zu erhalten, 

gekennzeichnet dadurch, dass 

f) die Verschiebungen der Phasenlagen aus Schritt c) registriert werden, 

g) Helligkeitsschwankungen der Beleuchtung erfasst werden 

h) Ausbleicherscheinungen des Objekts bei Fluoreszenzbeleuchtung bestimmt werden 
und die Verrechnung unter Beriicksichtigung der in den Schritten f) bis h) 
gewonnenen Ergebnisse erfolgt 

2. Verfahren zur Verbesserung der Tiefendiskriminierung optischer Systeme nach 
Anspruch 1, gekennzeichnet dadurch, dass Ausbleicherscheinungen des Objekts bei 
Fluoreszenzbeleuchtung dadurch bestimmt werden, dass aus vorzugsweise 
aufeinanderfolgend registrierten Helligkeitsverteilungen Gebiete gleicher 
Beleuchtungsintensitat bestimmt werden und die Helligkeit dieser Gebiete miteinander 
verglichen wird. 

3. Verfahren zur Verbesserung der Tiefendiskriminierung optischer Systeme nach 
Anspruch 2, gekennzeichnet dadurch, dass ein Normierungsfaktor durch eine 
Quotientenbildung aus den registrierten Helligkeiten der Gebiete gleicher 
Beleuchtungsintensitat erfolgt und mit Hilfe dieses Normierungsfaktor die 
registrierten Helligkeitsverteilungen normiert werden. 

4. Verfahren zur Verbesserung der Tiefendiskriminierung optischer Systeme nach 
Anspruch 1, gekennzeichnet dadurch, dass die Bestimmung der 
Objekthelliglceitsverteilimg m(x, y) durch die Gleichung 
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m ( x > y) = : — 

a 0 (x,y) 



erfolgt, wobei die Bestimmung des Vektors 



a = 



a x (x,y) 
\a 2 (x,y)) 



durch. L3sen eines Gleichuagssysterns aus der Verknupfung der registrierten 
HeUigkeitsverteilungen I j (x, y, mit den registrierten Phasenverschiebungen <j> i 
erfolgt, wobei i fur die Anzahl der Registrierungen der Helligkeitsverteilung steht 

5. Verfahren zur Verbesserung der Tiefendiskriminierung optischer Systeme nach 
Anspruch 3 und 4, gekennzeichnet dadurch, dass das Gleichungssystem die 
normierten Helligkeitsverteiluagen mit den registrierten Phasenverschiebungen 
verknilpft. 

6. Verfahren zur Verbesserung der Tiefendiskrimimerung optischer Systeme nach 
Anspruch 1, gekennzeichnet dadurch, dass eine Kalibrierung der Verschiebungen der 
Phasenlage erfolgt, wobei 

a) ein Objekt, welches vorzugsweise eine Spiegelflache ist, mit einer periodischen 
Struktur beleuchtet wird, 

b) eine erste Helligkeitsverteilung registriert wird, 

c) die Phasenlage der periodischen Struktur urn einen bestimmten geringen Betrag 
verandert wird, welcher registriert wird , 

d) eine zweite Helligkeitsverteilung registriert wird, 

e) der Unterschied der beiden Helligkeitsverteilungen bestimmt und bzgl 
streifenformiger Verteilung bewertet wird, 

f) die Schritte c) bis e) solange wiederholt werden, bis die Bewertung des 
Unterschiedes zwischen erster und folgender Helligkeitsverteilung auf ein Extremum 
fuhrt, und 

g) der so geflindenen Wert der Phasenlage der periodischen Struktur registriert wird. 
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7. Verfahren zur Verbesserung der Tiefendislmminierung optischer Systeme nach 
Anspruch 1, gekennzeichnet dadurch, dass durch nichtsinusformige Verteilung der 
Beleuchtung des Objekts auftretende Fehler korrigiert werden. 

8. Verfahren zur Verbesserung der Tiefendiskriminierung optischer Systeme nach 
Anspruch 7, gekennzeichnet dadurch, dass mittels Bandpassfilter hohere Harmonische 
der Grundfrequenz der periodischen Struktur aus der Helligkeitsverteilung 
herausgefiltert werden. 

9. Verfahren zur Verbesserung der Tiefendiskriminierung optischer Systeme nach 
Anspruch 1, gekennzeichnet dadurch, dass die Verrechnung der registrierten 
Helligkeitsverteilungen unter Beriicksichtigung von Korrekturwerten vorzugsweise fur 
die Beriicksichtigung der HeUigkeitsschwankungen der Beleuchtung, des 
Ausbleichens des Objekts bei Fluoreszenzbeleuchtung und nichtsinusfSrmiger 
Verteilung der Beleuchtung des Objekts erfolgt, welche durch lineare Optimierang 
bestimmt werden. 

10. Verfahren zur Verbesserung der Tiefendislaiminierung optischer Systeme nach 
Anspruch 9, gekennzeichnet dadurch, dass die lineare Optimierung auf eine Merit- 
Funktion der Form 

M(0 r9 d;b)^a o \F{a} o \ 2 +a x \F{a^ 

angewandt wird, wbbei F{5} eine Fuiiktionaltransformierte des Vektors 5ist, 0 t 
skalare Faktoren zur Charakterisierung der Beleuchtungshelligkeitsschwankung ist, 
d ein MaB fur dasr Ausbleichen des Objekts bei Fluoreszenzbeleuchtung ist, b ein 
Faktor zur Charakterisierung der niehtsinusf5rmigen Verteilung der Beleuchtung des 
Objekts ist und die a t Gewichtekoeffizienten zur Anpassung an die 

Aufoahmebedingungen vorzugsweise bzgl. verschiedener Signal-Rausch Abstande 
oder bevorzugter Frequenzen sind. 
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Zusamrnenfassung 



Die Erfindung betrifft ein Verfahren zur Tiefendiskrimiiiierung optisch abbildender 
Systeme. Es ist insbesondere in der Lichtmikroskopie anwendbar um erne Erhofrung der 
Bildqualitat bei der Untersuchung dreidimensional ausgedehnter Objekte zu erreichen. 
Sie ist anwendbar beim Verfahren der strukturierten Beleuchtung wie es in WO 97/6509 
beschrieben ist. Dazu werden Einfliisse durch Helligkeitsschwankungen der Lichtquelle, 
der Positionierung der abgebildeten periodischen Struktur und Ausbleichen des Objekts 
bei Fluoreszenzbeleuchtung ermittelt und bei der Berechnung der Objektstruktur 
beriicksichtigt 



Fig. 1 
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